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ABSTRACT 

The spectrum of radiation emitted following shock breakout from a star's surface with a power-law den- 
sity profile p oc x" is investigated. Assuming planar geometry, local Compton equilibrium and bremsstrahlung 
emission as the dominant photon production mechanism, numerical solutions are obtained for the photon num- 
ber density and temperature profiles as a function of time, for hydrogen-helium envelopes. The temperature 
solutions are determined by the breakout shock velocity vo and the pre-shock breakout density po, and depend 
weakly on the value of n. Fitting formulas for the peak surface temperature at breakout as a function of vo and 
po are provided, with 7^ ea k ~ 9.44exp[12.63(vo/c) 1//2 ] eV, and the time dependence of the surface temperature 
is tabulated. The time integrated emitted spectrum is a robust prediction of the model, determined by 7^>eak and 
Vo alone and insensitive to details of light travel time or slight deviations from spherical symmetry. Adopting 
commonly assumed progenitor parameters, breakout luminosities of sw 10 45 erg s _1 and ra 10 44 erg s _1 in the 
0.3-10 keV band are expected for BSG and RSG/He-WR progenitors respectively (7^ e ak is well below the band 
for RSGs, unless their radius is ~ 10 13 cm). > 30 detections of SN1987A-like (BSG) breakouts are expected 
over the lifetime of ROSAT and XMM-Newton. An absence of such detections would imply that either the typ- 
ical parameters assumed for BSG progenitors are grossly incorrect or that their envelopes are not hydrostatic. 
The observed spectrum and duration of XRF 080109/SN2008D are in tension with a non-relativistic breakout 
from a stellar surface interpretation. 
Subject headings: radiation hydrodynamics — shock waves — supernovae: general 



1. INTRODUCTION 

In a supernova (SN) explosion, a radiation mediated shock 
(RMS) traverses the stellar envelope. When the shock breaks 
out from the surface of the star, radiation is emitted with a 
prompt bright flash. Observations of breakouts provide im- 
portant information about the progenitor. A detailed under- 
standing of the expected bolometric and spectral properties 
is required for an accurate estimate of the progenitor's radius 
and for planning of transient search strategies. 

The bolometric and hydrodynamical properties of a non- 
relativistic planar R MS breaking out fro m a stellar envelope 
have been studied in ISapir et all d2011l) (hereafter Paper I). 
The analysis was preformed by neglecting matter thermal en- 
ergy and approximating radiation energy transport as diffu- 
sion with constant opacity k. The bolometric light curves and 
the hydrodynamic profiles were solved and fitted by analyti- 
cal formulas. T he observed proper ties of SN shock breakouts 
were derived in Kat z et all (120121) (hereafter Paper II) using 
the planar results, assuming spherical symmetry and taking 
into account limb darkening. 

The spectral energy distribution of the emission is depen- 
dent on the interaction of radiation and matter. In local ther- 
mal equilibrium (LTE), the radiation temperature is directly 
determined from the radiation energy density, and thus from 
the hydrodynamic solution. However, in general the radia- 
tion spectrum is determined by processes of absorption, emis- 
sion, scattering and adiabatic compression. As these pro- 
cesses have different time scales, the radiation state is time- 
dependent, and it may not be in LTE during shock propaga- 
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tion in the envelope. For instance, in fast RMSs there is not 
enough time for photons to be produced or absorbed in the 
shock front, and the radiation will reach LT E only in the far 
downstream JWeaverill976tlKatz et alj|2010l) . 

Earlier work studied the spectral properties of radiation 
emitted following breakout, either with numerical calcu- 
lations of specific progenitor s t ructure s (Klein & Chevalier 
1978; Ensman & Burrowsll 19921: iBlinnikov et al.lll 998. 2000; 
Tolstov 2010) or by analytical estimates (Katz et al. 2010; 
Nakar&Sari 2010). Past numerical work, however, either 
assumed LTE or used single-group approximations or multi- 
group calculations without fully accounting for the effect of 
Compton scattering on the spectral energy distribution. In 
particular, these works ignored incoherent Compton Scatter- 
ing, in which a photon changes its energy upon scattering 
on an electron. This effect plays a key role in the prob- 
lem, and drives radiation and matter into a Compton equi- 
librium. Consequently, the radiation temperature can be or- 
ders of magnitude higher compared to the LTE case (Weaver 
1976; Kat z et alll2010l) . Past analytical work applied either 
steady state solutions to the time dependent problem or as- 
sumed only local production of photons, not accounting for 
the role of diffusion in determining the emitted spectrum. 

This work follows the effective photon approximation used 
in the temperature calculation of steady state RMSs (Weaver 
1976; Katz et al. 2010), and applies it to solve the time de- 
pendent problem of planar shock breakout from a stellar en- 
velope. The problem is solved under the assumptions of lo- 
cal Compton equilibrium and bremsstrahlung as the dominant 
emission rate. Solutions are provided for the temperature and 
the emitted radiation spectrum. 

This paper is organized as follows. The formulation of the 
problem and the governing equations are presented in §|2] Re- 
sults for the time dependent profiles of the temperature and for 



the surface temperature are described in §|3] The surface tem- 
perature is then used to estimate the emitted spectrum in § [4] 
The expected emission from possible progenitors and the im- 
plied detection rates are discussed in § Comparisons with 
previous results in the literature and with the observations of 
the X-ray transient 080109/SN2008D are given in § The 
main conclusions are summarized in § [7] 

2. FORMULATION OF THE PROBLEM 

Consider the problem of a planar non relativistic RMS 
propagating in a stellar envelope, and breaking out from its 
surface. Assuming the energy density is dominated by radi- 
ation, and the opacity is dominated by Thomson scattering, 
the hydrodynamic problem can be described with flow con- 
servation equations coupled to a diffusion equation with con- 
stant opacity, and was solved in Paper I, In order to solve for 
the radiation temperature, further assumptions are required: 
here we solve for the temperature assuming the radiation is 
in local Compton equilibrium and the emission mechanism is 
bremsstrahlung emission, 

The initial density profile is assumed to be a power-law of 
the distance from the surface, p oc x". In terms of the opti- 
cal depth with respect to the surface, r = nfpdx, the density 
profile is given by 



p(r) = po(/3,,rr /( " +I) , 



(1) 



where /3o = vo/c. At large optical depth, the asymptotic shock 
velocity can be described by the pure hydrodynamic solution 
(Sakurai 1960), and parameterized by 



v s h- 
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(2) 



The values of po and vo are the density and velocity at the 
point r = c/v s h, which is obtained in a pure hydrodynamic 
shock propagation. 

In Compton equilibrium, the radiation is characterized by 
a Wien spectrum, dejdv oc i/ 3 e~ v / T ~<, with radiation temper- 
ature, T 7 , given by the electron temperature, !T 7 = T. The 
relation between radiation temperature, pressure and photon 
number density is 



n 1 



(3) 



Note that in the case of LTE T w 0.9p/« 7 , and eq. (0 intro- 
duces a 10% error in the radiation temperature. 

The number density of photons is determined by photon 
transport in space and photon production by the plasma, 
which depends on the plasma temperature and density. As the 
hydrodynamic flow equations are independent of temperature, 
the detailed interaction of photons with matter does not affect 
the hydrodynamics. The independent solution for the energy 
density and density profiles, found in Paper I, can therefore be 
used to calculate the photon number density and temperature 
profiles as a function of time. 

2. 1 . Hydrodynamics and energy diffusion equations 

We briefly describe the hydrodynamic equations which are 
given in Paper I. Lagrangian coordinates are used to describe 
the flow. The mass coordinate m equals the mass per unit area 
enclosed between each mass element and the surface. The 
equations describing the position, x, velocity v and energy 



density e as a function of m and time t are 

d t x = v, 

d,v = -d m p, 

d,(e/p) = -d m j e -pd m v, 



(4) 
(5) 
(6) 



where p = d x m is the density, j e = -c(3K) _1 9 m e is the energy 
flux (energy current density), p = e/3 is the radiation pressure 
and k is the constant opacity. The boundary conditions at the 
surface are taken as 

e (/3 r = 0) = 0, d t v(PoT = 0)=~j e (/3 T = 0). (7) 

c 

The dimensional parameters n,po and vo along with the di- 
mensionless parameter n completely define the hydrodynamic 
problem. The hydrodynamic profiles can be obtained, solv- 
ing the equations of motion ©-(O, using the initial density 
profile (U), the asymptotic shock velocity (f2]i and the surface 
boundary conditions dT) (Paper I). 

2.2. Photon number density and temperature 

In the effective photon approximation the photon number 
density is calcul ated from the photon emission rate, corrected 
for absorption (Weaver 1976). Although the photon emission 
rate diverges logarithmically at low energies, photons that do 
not up-scatter in energy are absorbed, and a lower cutoff en- 
ergy can be introduced. This cutoff energy is the lowest en- 
ergy from which a photon can be "thermalized" to the plasma 
temperature T. 

The transport of photons is modeled by a diffusion equation 
with an additional source term, Q 1 , which accounts for the 
local net production rate of photons, 

d,(n 7 /p) = -d„,j 1 + Q 1 , (8) 

where the photon number current is 

c 

j~ f =-—dmi~ f . (9) 

Note that Compton scattering conserves the number of pho- 
tons, and therefore scattering does not directly appear in the 
equation for the photon number density. 

At large optical depth the radiation is found at LTE. The 
inner boundary condition for the photon number current is 
accordingly taken as 



7 7 (/? T ->• oc) =- 



e(/3 r -> oo) 
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M(3 T -»■ 00), (10) 



where abb = (tic)~ 3 ir 2 /l5 is the Stefan-Boltzmann energy den- 
sity coefficient and the 0.9 factor accounts for the error intro- 
duced using eq. (01 in LTE. At the surface, the photon number 
density boundary condition is taken similar to the energy den- 
sity boundary condition, 



« 7 (A)T = o) = o. 



(ii) 



Note that the boundary conditions at the surface are taken as 
zero energy and zero photon number density (see also Paper 
I), but the photon temperature is finite. 

The source term in eq. (0 is the net photon production rate 
per unit mass, which depends both on temperature and density 
of the plasma, 



1 

— ) 

P 



Q~ l (pJ)= Z qf(pJ)UUp-T), 



(12) 



where q?r is the photon number density emitted per unit time 
and / a b s is a dimensionless factor which accounts for absorp- 
tion. 

In the temperature range 10 eV < T < 50 keV, the dominant 
mechanism for photon production in a fully ionized hydrogen- 
helium plasma is electron-ion bremsstrahlung emission, while 
pair production and bound-free emission are negligible. Con- 
sider a completely ionized plasma with atomic number Z = 
E,jc,'Z,' and atomic weight A = £,\x,-A/, where x,- are the atomic 
fractions of each element in the envelope composition. The 
effective photon emission rate is given by dWeaverj 119761; 
iKatzetal 



photon 
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(13) 



where ctt is the Thomson cross section and the factor A ge ff ex- 
presses the integral of the bremsstrahlung emission spectrum 
down from a low energy cutoff v c (as this integral diverges 
logarithmically). The integral over the emission spectrum, in- 
cluding the gaunt factor calculated in the Born approxima- 
tion, can be described as a product of t wo separate fu nctions, 
Ag.eff = gi(A)£i(A), where A = v c /T (IWeaveri 119761) . These 
functions can be calculated numerically. The function g\(X) 
is approximated by 

gi(A)= 1.226-0.4751og(A) + 0.00131og 2 (A), (14) 

and E\{\) is the exponential integral function 



£l(A): 



-dx, 



(15) 



which can be approximated as £i(A) = -log(A)- 0.5772 for 
small A. 

The cutoff energy v z is determined by the lowest energy of 
emitted photons that reach thermalization with the plasma by 
up-scattering, before being absorbed. The energy of photons 
which up -scatter T/(m e c 2 ) times before getting a bsorbed is 
given by (Cha pline & Stevenslll973tlWeaveril 19761) . 
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(16) 
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where p = 10 9 p-9 S cm 3 an d T = lT^eV keV. When the 
Compton y-parameter is lower than unity, T /(m e c 2 )T 2 < 1, 
photons escape the medium before being up-scattered in en- 
ergy. Therefore, A g e ff has a minimal value of unity. 

In order to account for absorption, t he following te rm is 
introduced in the effective emission rate (Weaver 1976), 



/abs = 1 ' 



(17) 



where the equilibrium temperature is calculated from the local 
energy density, 



T = 



0.9flP 



1/4 



(18) 



This ensures that there is no net photon production when the 
radiation is in LTE with matter. 

To conclude, the photon temperature and number density 
as a function of time are described by the photon diffusion 
equation ^ with the boundary conditions (TToT > and (fTTT i. and 
using the Compton equilibrium relation (01, the effective pho- 
ton generation rate dl2t and the hydrodynamic profiles that 
are obtained independently (see § 12. 11 1. 

2.3. Dimensionless parameters 

In the considered temperature range, 10 eV < T < 50 keV, 
and density range, p < 10~ 6 g c m" 3 , opacity is dominated 
by electron scattering (ICoxl fl968). The total opacity is then 
k = (£,'X,Z;/£,jic,'A,')(ot/'Wp)- Comparing the radiation energy 
density emitted in a dynamical time to the shock energy den- 
sity, q^Ttd = pv'g, the following photon temperature scale is 
obtained 



To = 



^iXfAi 



Plm v c 2 



Introducing dimensionless variables, a = a/ao, with 



(19) 



k) = 5, e = pov , n 7 .o = eo/7b, m = (k/3q) , (20) 

k-PoVq 



eq. dD can be re-written as 



e Ti 



djn 7 = —dfp+-d lT ,n 1 + qt / I 1- - 4 jA 



(21) 



with a dimensionless emission rate 
q e J l (p,T)=^ T77 a 



2 3 / 2 fmX l2 ( pVfT^- 111 
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(22) 



(23) 



Note that i/ c , and therefore also A ge ff, scales with 

For non-relativistic velocities, a single dimensionless pa- 
rameter appears in the hydrodynamic problem, namely, the 
power-law index n. Thus, for a given index n, the hydrody- 
namics follow a unique solution which is scalable with the 
parameters K,po and vo (see Paper I). However, the tempera- 
ture solution does not generally produce such a universal be- 
havior, as there are two more dimensionless parameters in the 
problem. 

Eq. (l2Tb implies the existence of three interesting regimes. 
At LTE the temperature behavior is directly determined by 
the universal hydrodynamic solution, and therefore T/T eq fi 
exhibits a universal behavior. When v c > 0.45 T, A ge ff is at 
its lower limit of unity, and there are two energy scales in the 
problem in addition to the shock energy. The first is the equi- 
librium temperature r eq ,o and the second is the temperature 
obtained in bremsstrahlung emission over a dynamical time, 
Tq. For a given index n, the photon number density and tem- 
perature in this regime depend on the dimensionless parame- 
ter To/T eq o <x pn 0o ■ When T 3> T eq o there are again two 
energy scales in the problem (in addition to the shock energy), 
Tq and the lower cutoff energy v c , and the solution is deter- 
mined by 2 dimensionless parameters, n and T /v c ex p /3q S . 



3. RESULTS 

The photon diffusion equation, eq. ([8j, was numerically 
solved using an implicit scheme in time, on a uniform mass 
grid. The hydrodynamic equations, eqs. dUi-© were solved 
simultaneously, in a manner described in Paper I, In order to 
achieve fast convergence, the initial photon number density 
was derived from the initial energy density assuming LTE, 
n-y = (0.9aBB) 1 ' 4 e 3 ' 4 /3. The inner boundary condition, eq. 
( flOb . was obtain ed from Sakura i's time dependent hydrody- 
namical solution (Sakurai 196U- 

The resulting profiles of the temperature as a function of op- 
tical depth at different times are given in § 13.1 I for a hydrogen- 
helium envelope with n = 3. The ti me d ependence of the sur- 
face temperature is described in § 13.21 Note that as calcu- 
lated temperatures scale with atomic number and weight as 
S,x,Z?/S;x,A, (see eq. (fT9b ). the following results are appli- 
cable to envelopes with any hydrogen-helium composition. 

3.1. Temperature profiles 

Temperature profiles at different times during shock prop- 
agation, calculated with an ion number density, n Pj o = 
po/{Am p ), of fUfl = 10 I5 A _I cm" 3 and different values of /3o 
are shown in figureQ] Also presented are the black-body tem- 
perature profiles, calculated from the energy density profiles 
using eq. (fl8t , As is presented, temperatures in excess of 
2 keV are obtained for /3o = 0.2. In addition, even for a low 
shock velocity, (3q = 0.05, the temperature at an optical depth 
of (3q 1 is a factor 2 higher than the LTE value. 

There are several competing effects determining the tem- 
perature at the shock front, which is affected both by pho- 
tons produced locally at the shock front and by photons dif- 
fusing from the downstream. Diffusion of photons from the 
far downstream lowers the radiation temperature, as there are 
more photons to share the shock's energy. The number of pho- 
tons required to achieve thermal equilibrium decreases as the 
shock propagates outwards and its energy density decreases. 
On the other hand, bremsstrahlung photon emission decreases 
as the shock propagates outwards to lower densities, and the 
optical depth of the post-shock shell from which photons may 
diffuse up to the shock, At = 7/3~ h ', is decreasing as the shock 
accelerates. 

LTE is obtained at the shock front when there is enough 
time for photons to be produced at the front, or reach it by dif- 
fusion from the downstream. In the far downstream the radia- 
tion is at LTE, but as the shock traverses the gas, the radiation 
is driven out of thermal equilibrium at the shock front. This 
occurs also for low shock velocities, in which LTE is naively 
expected to hold even close to breakout. In a steady state 
shock, mass elements in the downstream move with the same 
velocity, and the shock front is affected by photons diffusing 
from an optical depth of 7/3q' into the downstream. However, 
as the shock gets closer to the star's surface the matter velocity 
gradient steepens, and matter at the shock front moves faster 
than matter in the downstream. As a result, the shock front 
"sees" a smaller optical depth into the downstream than in the 
steady-state case. 

It is possible to roughly estimate the range of /?o and 
Pa for which L TE is achieved from the steady-state RMS. 
IWeaverl (11976b concludes that the shock front of a steady- 
state RMS would be in LTE for shock velocities (3q < 
0.02(po/10~ 9 g cm" 3 ) 1 / 30 , based on the effective photon ap- 
proximation also employed in this work. Nakar & Sari (2010) 



conclude that the shock front would be in LTE for shock ve- 
locities (3q < 0.05(po/10~ 9 g cm" 3 ) 1 ' 30 , based on local pro- 
duction of photons in the immediate downstream. O ur results 
show that the estimate provided by IWeaverl (119761) is more 
accurate for the shock breakout problem. 

Analytical estimates for the maximal temperature at the 
shock front can be obtained by applying the steady-state RMS 
temperature solution to Sakurai's hydrodynamical solution. 
Taking the initial density and shock velocity values as a func- 
tion of optical depth from eqs. (H)-I0, and using them as 
a uniform input for an RMS propagating in a homogeneous 
medium (IWeave r 1976; Katz et al] [2010). allows one to esti- 
mate the maximal temperature as a function of optical depth. 
A comparison with the steady-state estimate reveals good 
agreement at large optical depth, but an increasing discrep- 
ancy as the shock approaches the surface. At an optical depth 
of r = /Sy 1 the calculated temperature may a factor 2-5 lower 
than the analytical estimate. This difference originates from 
the large photon number density in the downstream, compared 
to the steady-state case. Thus, more photons arrive from the 
downstream, lowering the temperature at the shock front. As 
the shock accelerates and radiation is driven away from equi- 
librium, the effect becomes more pronounced. Hence, cal- 
culating the emitted spectrum from the steady-state estimate 
may produce large errors. 

Figure|2]presents temperature profiles at different times fol- 
lowing breakout, calculated with different values of /3rj and po- 
Photons arriving from an optical depth r = /3q scatter on their 
way to the surface, but almost no new photons are produced. 
As n 1 >■ Zn v , the electrons' temperature is expected to be 
equilibrated with the photons. Thus, the surface temperature 
following breakout is similar to the temperature at r = /3q : . 
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FIG. 1. — Temperature profiles as a function of optical depth at different 
times prior to breakout (solid lines). Also plotted are the LTE temperature 
profiles (dashed line). fSgc is the shock velocity "at breakout", i.e. the velocity 
for which t = /3r' in the Sakurai solution, and tg is defined in eq. j20> . 



3.2. Surface temperature 

The observed emission is determined by photons escap- 
ing from the surface of last scattering, which is at an optical 
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FIG. 2. — Temperature profiles as a function of optical depth at different 
times following breakout (soli d lines). Also plotted are the black-body tem- 
perature profiles from eq. ( 1181 (dashed line). 

depth of t w 1 . The emitted radiation spectral energy distri- 
bution can be described using the surface temperature, T(t) = 
T(Pqt — > 0,f), in the limit /?o *C 1. When Thomson opac- 
ity dominates over bremsstrahlung opacity, radiation is not at 
LTE at this optical depth and the surface temperature does not 
equal the effective temperature, T e s = [4£/(aBBc)]'^ 4 , where 
C is the bolometric emitted flux. Thus, for the scattering- 
dominated case the surface temperature is very close to the 
temperature at an optical depth of r = /3q' (see figure|2]i. This 
is due to the fact that photons with energies < [3Q 2 m e c 2 do not 
lose significant energy diffusing from r = (3q 1 to t = 1, and 
that the energy required to raise the electrons' temperature to 
the radiation temperature is small. 

The surface temperature rises sharply, with a peak width of 
~ to, and then drops with time. The peak surface temperature 
7^eak as a function of the breakout parameters /3q and « p o is 
shown in figure|3]for n = 3 and n = 3/2. Defining (3q = 0.1/3_i, 
the results can be fitted by the formulas 



+ (o.26-0.08/3±{ 2 ) log 10 (A« p ,i 5 ) , n = 3 



logi ° \^T) 



log 



10 



/peak 

~eV~ 



•69/3lf 



=0.95 + 1.78/3l{ 2 



+ (0.26-0.08/3X) log 10 (An pA5 ) , n = 3/2 

(24) 

with an error of less than 10% in 7^ e ak in the range 0.4 > /3n > 
0.01 and 10 17 A-' cm" 3 > n pfi > l0 n A~ l cm" 3 . The results of 
the fitting formulas averaged between n = 3 and n = 3/2 are 
also shown in figure [3] in dashed lines. The averaged fitting 
formula is given by 



ql/2 



log 1( ,(^) =0-975 +1.735/3l{ : 

+ (0.26-0.08/?l{ 2 )log 1() (A» p , 15 ). 



(25) 



The peak temperature time is roughly f pea k = — Ito 
prior to the peak luminosity time, which is f re f/fo = 
-1.25rn (-0.78f ) for « = 3 (« = 3/2) (see Paper I). The peak 
time can be described as a function of the parameter \ = 
/3 (An P0 /10 15 cm" 3 ) -1 / 36 and fitted by the formulas 



fpeak =fref" (0.83 - 0.85e" 141x ) t Q , R = 3 
fpeak =?ref- (0.7 - 0.8e" 1775x ) f , B = 3/2 



(26) 



with an error of less than 10% in f pea k in the range 0.3 > x > 
0.05. 

Figure |4] shows the surface temperature normalized to the 
peak value as a function of normalized time for several cal- 
culations with different values of vo and « Po = 10 15 A _1 cm" 3 . 
The asymptotic behavior of the surface temperature is shown 
to be bounded between a power-law slope of oc t~ ] /3 and a 
power-law slope of oc f " 2 / 3 . The temperature curves are also 
provided in table[TJ(|2]l for n = 3 (n = 3/2). The derived spectral 
luminosity is a robust prediction of the shock breakout model, 
since both the temperature and the bolometric luminosity are 
not sensitive to the envelope density structure. 

A T oc f" 1 / 3 asymptotic decline is expected for radiation in 
thermal equilibrium or in adiabatic expansion, in agreement 
with the power-law upper bound presented in figure 2] A 
T oc f" 2 / 3 asymptotic decline was predicted by iNakar & Saril 
(2010) for the out-of-equilibrium case, in agreement with the 
exact solution presented here. We note however that at late 
times, the production of additional photons turns out to be 
insignificant as opposed to the assumptions leading to the 
derivation of T oc r 2 / 3 in INakar & Sanl (120101) . 

Different combinations of vo and pn, which give the same 
value of x, yield approximately the same temporal depen- 
dence of the surface temperature (normalized to T pei k, see 
§ 12.31 for discussion). The time dependence of the surface 
temperature is given in tables Q] and |2] for various values of 
X- In order to calculate the normalized surface temperature at 
different times for values of x not provided in the tables, an 
interpolation in x is required. A linear interpolation, 

T(t,x) T(t, X i) , fT(t, X 2) T(t,xi)\ x-xi 



7peak(x) 7peak(Xl) \7peak(X2) 7^,eak(Xl)/ X2~Xl 

(27) 
where xi > Xi are adjacent values provided in the tables and 
X2 > X > Xi , yields an error of less than 5% in the normalized 
temperature. Eq. (l24l can then be used to calculate the surface 
temperature in physical units. The time in physical units is 
obtained by using eq. (1201 ) and the peak surface temperature 
time which is approximately given by eq. d26l l. 

4. EMISSION SPECTRUM 

Photons diffusing outwards from an optical depth of r = /3q ' 
are found in local Compton equilibrium and are characterized 
by a Wien spectrum close to the peak energy. The radiation 
spectrum is not significantly modified while photons diffuse 
to the surface, as the electron temperature is equilibrated with 
the radiation. Therefore, the escaping photons' spectrum can 
be estimated from the bolometric emitted flux and the surface 
temperature, assuming a Wien spectrum. Note that no high 
energy power law tail is expected in the radiation spectrum, 
as the radiation temperature is higher than or similar to the 
matter temperature. 

The instantaneous emitted energy flux per logarithmic unit 
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photon energy is given by 



C u (t) = -C(t)x exp(-x), x = v/T(t) 



(28) 



where the pre-factor 1 /6 is the normalization factor, so that 
J_ C v {t)dh\v = C(t). The instantaneous emitted spectral flux 
is used t o der ive the spectral fluence (time integrated spectral 
flux) in § 14.11 and the observed luminosit y in d ifferent energy 
bands for spherically symmetric SNe in § [ 



4. 1 . Fluence spectrum 

While the observed spectrum of a SN shock breakout may 
be sensitive to light travel time averaging and to deviations 
from instantaneous arrival of the shock to the surface at all 
angles, a robust estimate can be obtained by considering the 
fluence spectrum: the flux spectrum integrated over time, 



£#)= / Cv(t')dt'. 



(29) 



The flux spectrum integrated up to time (t— ?ref)Ao =100 for 
different values of po and vo are shown in figure \5\ Normaliz- 
ing the emitted fluence by 



£o = k r A)C 2 



(30) 



and the photon energies by the peak surface temperature, an 
almost universal shape is obtained for the fluence spectrum. 
Therefore, the photon energy of the integrated flux peak is a 
good estimate for the maximal surface temperature, 



^peak — 3 'peak- 



(31) 



A relation between shock velocity and density at breakout can 
then be inferred from the observed fluence peak energy and 
eqs. (|24) and ( T3Tb . 

The low energy part of the integrated spectrum can also 
be estimated from eqs. (|28| | and ( l29l . Taking T oc t~ s and 
C oc f" 4 / 3 , the low energy tail of the integrated spectrum 
has a non-thermal part with a photon energy dependence of 
Sy oc i/V(3*), This result can also be derived by c onsidering 
e v oc Cj and taking T ~v dNakar& Saril 1201 Oh . Thermal 
(Rayleigh-Jeans) contribution to the low energy tail of the 
instantaneous flux spectrum, with a v^ dependence, would 
not change the integr ated spectrum significantly, assuming 
1/3 < S < 2/3 (see § 13.21 ). although a more moderate spec- 
tral slope may change this result. 

We point out that the fluence spectrum is different from the 
instantaneous emitted flux (Wien) spectrum as well as from 
a black-body spectrum. In particular, while a Wien spectrum 
peaks at ^ pea k = 47^ e ak and a black-body spectrum peaks at 
t'peak = 2.827peak, the fluence peaks at i/ pe ak = 37~, e ak- In addi- 
tion, the low energy tail at energies v < 3T ped k has a shallower 
slope than either Wien or black-body spectra. 

4.2. Luminosity in different energy bands 

Considering shock breakout from a spherically symmetric 
SN and taking into account light travel time integration, the 
observed flux in different energy bands can be derived follow- 
ing the derivation of Paper II for the bolometric luminosity, 
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(32) 
where R is the stellar radius, hip) is the angular distribution of 
radiation intensity at the surface, V\ and V2 are the energy band 
boundaries and F a bs(Ti/,abs) is tne absorption fraction, where 
Tj/,a is me absorption optical depth on the line of sight as a 
function of photon energy. 

Assuming that the radiation angular distribution is not fre- 
quency dependent and ignoring absorption, the integral over 
photon energies can be solved separately, 
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ton energy, from calculations with different values of breakout shock velocity 
/?o, ion number density n p o = 10 15 A~' cm -3 and n = 3. 

where t' = t-R(l-fi)/c, and the integral T = J Xl x 4 exp(-x)dx 
is given by 

JvOi ,V2,T)= e~ x (x 2 + 3x 2 + 6x+ 6) |^'/r ■ < 34 > 

Therefore, whenever the surface temperature crosses in- 
side/outside an energy band boundary, the luminosity in this 
band lights up/drops exponentially. 

5. SHOCK BREAKOUT SEARCHES 

Assuming all SNe are accompanied by a shock breakout, 
detecting the shock breakout transients and their signature can 
be used to reveal the properties of SN progenitors. In order to 
plan a search strategy for shock breakout transients it is nec- 
essary to estimate the expected luminosity and photons' en- 
ergy. Based on the exact results presented here and in Papers 
I- II, th ese values are estimated for possible SN progenitors in 
§ I5.ll followed by estimates for the detection rate o f sho ck 
breakout events for a number of X-ray telescopes in § I5.2I 

5.1. Possible progenitors 

We focus next on red supergiants (RSGs), blue supergiants 
(BSGs) and He-rich Wolf-Rayet (WR) stars as possible pro- 
genitors of core collapse SNe. RSGs and He-rich WRs have 
been associated wit h SNe of type II-p and type lb respectively 
dSmartt et al.l 120091) . while the progenitor of SN1987A has 
been identified as a BSG with a mass of ~ 20M© and radius 
of ~ 3 x 10 12 cm dArnett et al.lll989l) . 

From pre-SN images, the mass range of progenitors for 
type II-p SNe is constrained to ~ 10-20M Q , while the pro- 
genitors' radii are i?~3-7xl0 13 cm dSmartt et all l2004t 
iMaund et al.ll200l: iMattila et al J 120081 iVanDvket alJl2012h . 



Using an analytical model to infer the radius of SNellp 
progenitors, different values were obtained by different au- 
thors, R«0.5-4x 10 13 cm by iHamuvl ( 120031) and R w 
1.4-9 x 10 13 cm by iNadvozhir] ( 120031) . Both ranges are 
within the known radii distribution of RSGs (Ivan Belle et al.l 
2009). Compared to the supergiants, WR are compact with 



R ~ 10 1 ' cm and have lost their hydrogen shell (lLangerll989l) . 
In particular, helium rich WRs have been suggested as pro- 
genitors of type lb SNe. 

The physical parameters of the shock breakout problem, 
the breakout shock velocity and pre-shock density, can be es- 
timated from the stellar parameters of each progenitor type, 
yielding (see equations 14 and 15 in Paper II) 



p = Q.Q5M% l3 vlf 5 K 
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and 
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(36) 



where the envelope mass is M e j = xM x Mq, the star radius is 
R = KYRx cm, the typical ejecta velocity is v* = 10*v* , x km s _1 , 
and f p is a dimensionless parameter (see Paper II for discus- 
sion). Note that the trend in shock velocity of the different 
progenitors is consistent with the trend in the photospheric 
velocities of the different SN types, SNellp < SN1987A < 
SNelb (evaluated from observations made late in the spheri- 
cal stage, after hydrogen/helium recombination). 

Assuming that the stellar envelopes are composed of ion- 
ized hydrogen and helium, with mass fractions X and Y re- 
spectively, the opacity for each progenitor is taken as k = 
(l+X)/5 cm 2 g _1 . For BSGs and RSGs (WRs), these mass 
fractions are taken as X = 0.7 and Y = 0.3 (X = and Y = 1). 
Using eqs. (A1)-(A2) in lCalzavara & Matznerl (120041) and as- 
suming that the core mass is small compared to the envelope 
mass, the parameter f p can be calculated for BSGs and WRs, 
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(37) 



where // = (2X + 0.75F) -1 is the mean molecular weight and 
L* = 10 A L V L© is the stellar luminosity. For RSGs the pa- 
rameter fp is not well-constrained , and we take its value to 
be unity (see also discus sions in Matzner & McKee 1 19991 : 
ICalzavara & Matznerj|2004l) . 

In addition, the typical ejecta velocity can be estimated as 
v* = CEjn/Mej) 1 / 2 where E- m = WE X erg is the explosion en- 
ergy (Matzner & McKee 1999). Unfortunately, this parameter 
can not be inferred directly from observations, but only esti- 
mated to an order of magnitude. Given these relations, and 
using eqs. (124-b and OTI ). the shock breakout flash would peak 
at X-ray energies, f pea k > 0.5 keV, if the progenitor's radius 



is roughly (up to logarithmic corrections) 
R x <4.5 x 10 12 M7n' 67 £c 2 i 17 



1+X 
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cm (RSG) 
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1.7 
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I + Y^ 066 /2X+ 0.757 
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V 0.62 
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cm (BSG,WR). 

(38) 



As expressed by these equations, the maximal progenitor ra- 
dius for producing an X-ray flash at shock breakout strongly 
depends on both the explosion energy and the progenitor's 
mass. For instance, if SN1987A-like events are produced by 
BSGs with fi<5x 10 12 cm, it is expected that they would 
flash at X-ray energies. Additionally, RSGs with radius of 
R = 10 13 cm and explosion energy of E ln = 2 x 10 51 erg pro- 
duce an X-ray flash at shock breakout. While the progenitors' 
stellar parameter distribution is unknown for each SN type, 
this raises the possibility that some SNellp may be accompa- 
nied by an X-ray flash. 

For demonstration purposes, we consider the following stel- 
lar parameters: 



M ej = 10M Q , fi = 3xlO l3 cm 
M ej = 10M Q , R = 3x 10 12 cm, U 



10 5 Lr. 



(RSG) 
(BSG) 



M ej = lOM , R=l x 10 n cm, L* = 1O 53 L (WR), (39) 

and take for all progenitors an explosion energy of E m = 
10 51 erg. Estimating the peak breakout temperatures using 
eq. ( f24b yields 

7^eak = 49.4eV (RSG) 

7^eak = 258eV (BSG) 

T P eak = 3.7keV (WR). (40) 

For comparison, Naka r & Saril (120101) predict breakout 
temperatures of r RSG = 63 .5M^ 3 E^ 5 R^ 65 eV, T BSC = 



IM^ Z E^R-^ A keV and T WR = 25. SM^'E^fR-^ keV, 
yielding T = 3 1 . 1 eV, 301 eV, and 25.8 keV for the RSG, BSG 
and WR progenitors we are considering. Coincidentally, the 
peak temperature fo r the BSG prog enitor in our model and 
in the model by Na kar & Saril (120101) is similar, but the de- 
pendence of the peak temp e rature on energy is stronger in the 
estimates of Nakar & Saril (120101) than in our results. In the 
explosion energy range of £'i n = 0. 5-2 x 10 51 erg the spread in 
the peak temperature is a factor of 3 in our m odel and a factor 
of ~ 10 in the model by Nakar & Sari (2010). The differences 
between these results and ours lead to significant differences 
in the luminosities predicted in different bands and hence for 
outburst detectability, as discussed below. 

The luminosities in different energy bands following break- 
outs from RSGs, BSGs and WRs are presented in fig- 
ures |6][8] for the energy bands of the burst alert telescope 
(BAT) (15 keV < v < 150 keV), the X-ray telescope (XRT) 
(0.3 keV < v < 10 keV) and the optical/UV telescope 
(UVOT) (2 eV < v < 7 eV) on board Swift. The band lumi- 
nosities are calculated using eq. (l33l up to time f s = R/(4vq) 
when the planar approximation is no longer valid (see Paper 



II). For completeness, the b and luminosities in the sphe rical 
stage, calculated according to Rabinak & Waxmanl ( l201 ll) . are 
also presented in the figures. 

It is apparent from figures [7] and |8]that most of the prompt 
emission (on a timescale R/c) in BSGs and WRs is in the XRT 
band. Figure|8]also shows that for WRs there is also a smaller 
contribution to the emission by the BAT band. Note though 
that this result is sensitive to the assumed explosion energy, 
as higher values of this parameter imply that the emission 
would be dominated by the BAT band. In RSGs, although the 
emission following breakout peaks below the lower energy 
of the band, ^ pea k = 37J, ea k < 0.15 keV, there is still signifi- 
cant emission in the XRT band (see figure|6]l. This is a result 
of eq. (l34l . as the luminosity radiated in the XRT band is a 
fraction (x 3 + 3x 2 + 6x + 6)e _JC /6 of the bolometric luminosity, 
where x = v\ /T(t) and v\ is the lower energy of the band. 

Essentially, the total number of photons emitted in the burst 
determines the detectability of the transient. Approximating 
the emitted number of photons as iV 7 = 47r/? 2 £o/7^eak, yields 
for the different progenitors 



JV 7 =9xl0 3/ 
JV 7 =5.2xl0 55 
JV 7 = 1.6xl0 52 



(RSG) 
(BSG) 
(WR). (41) 



This approximation is good to better than 5% accuracy for 
an energy band between 7peak/3 <v< 57JSeak- Note that the 
larger the radius of the progenitor, the lower is the breakout 
shock velocity (see ( |35l l) but the higher is the energy released. 
Consequently, as the peak temperature is also lower, the to- 
tal number of emitted photons is higher, which makes the 
event easier to detect. However, as the peak temperature is 
lower, the emitted photons may be absorbed by the interstel- 
lar medium. Considering a neutral hydrogen column density 
of A^h = 1 x 10 21 cm" 2 , the optical depth for photoionization is 
> 1 for v < 0.5 keV dWilmset al. 2000), and photons emitted 
at lower energies are efficiently absorbed. 

Even though the emission following breakout from an RSG 
with radius R = 3 x 10 13 cm peaks at energy v ve! ± = 37peak < 
0.15 keV, the number of emitted photons is so high that it is 
still observable in the XRT band if N n < 1 x 10 21 cm -2 . This is 
also seen in figure|6] Using eqs. (1291 and ( T33l , the number of 
observed photons is diminished by a factor (x 2 + 2x+2)e~ x /2 
compared to the number emitted, where x = V\ /7peak and V\ 
is the lower energy of the band. However, in this case only 
the high energy tail of the fluence spectrum would be ob- 
servable, making a definite identification of the transient as 
a shock breakout event difficult. 

Thus, the best prospects for detecting a shock breakout 
flash are from BSG progenitors with radii R fa 3 x 10 12 cm, 
the burst from which is characterized by a relatively long 
duration, ~ R/c < 100 s, and peaks in low X-ray energies, 
*Vak = 37^eak ~ 0.8 keV. In addition, type II-p SNe may be 
accompanied by an observable X-ray flash if the neutral hy- 
drogen column density is Nn < 1 x 10 21 cm" 2 and the progen- 
itor's radius is R < 10 13 cm as proposed for some SNellp by 
lHamuvl(l2003l) . 

5.2. Detection of SN breakouts 

In general, detection of a transient is determined by the 
background and by the sensitivity of the detector. The X- 
ray background in the energy band 0.3 keV < v < 10 keV 



46 



10 



1(T 



H 10 

I 

M 

iio 40 



10 



38 



10 



36 























Bolometric 

BAT 

XRT 

UVOT 









10" 



icr io H 

t [s] 



10" 



FIG. 6. — Luminosity in different energy bands, for breakouts from an RSG 
associated with SNellp. The band luminosities in the planar stage (spher- 
ical stage, Rabinak & Waxman 1 20 1 1 )) are presented in solid lines (dashed 
lines). Colors (black, blue, red, green) represent the different energy bands 
(bolometric, BAT [15, 150] keV, XRT [0.3, 10] keV, UVOT [2,7] eV). 



46 



10 



10 



H 10 

60 
^ 40 

., 10 4 ° 



10 



38 



10 



36 




■ Bolometric 
BAT 
XRT 
UVOT 



10" 



10" 



t \s] 



FIG. 7. — Luminosity in different energy bands, for breakouts from a BSG 
associated with SN1987A. Notation is the same as figure|6] 



is »cxb ~ 001 counts cm 2 deg 2 s ' dDe Luca & Molendl 
12004 iHickox & Markevitchl 120061) . while in the energy 
band 15 keV < v < 10 keV it is n C X B - 7.5 x 
10" 4 counts cm" 2 deg" 2 s" 1 (lAjello et all 120081) . The num- 
ber of photon counts «d e t needed for a significant detection of 
a transient with duration Af with a detector with an effective 
area A e ff, PSF angular size $, field of view SI and lifetime t\ is 
determined by Poisson statistics, (a" del e~ a )/ndet = L , where 
a = »cxB(Aeff/cm 2 )($/deg 2 )(Af/s) and L = fo/Af)(ft/*). 
Because of the low X-ray background, «d e t < 10-25 for ac- 
tive and non-active X-ray telescopes (see table [3] for X-ray 
telescopes characteristics). However, «det only represents a 




FIG. 8. — Luminosity in different energy bands, for breakouts from a WR 
star associated with SNelb. Notation is the same as figure [6] 



lower limit on the number of photon counts needed for iden- 
tifying a transient over exposure of duration At. To identify 
the transient with SN shock breakout it is necessary also to 
compare the total emitted energy and the fluence spectrum of 
the event with the breakout model. We therefore arbitrarily set 
the lower limit on the number of detected photons at «d e t =100 
in order to determine the nature of the event according to the 
fluence spectrum. 

Assuming emission of /V 7 photons in a stated energy band, 
such events are detectable to a distance of 
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Taking the shock breakout rate per unit volume to be «b, the 
number of shock breakout events detected is 
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Table [3] presents the expected number of shock breakout de- 
tections over the lifetime of several active and non-active X- 
ray telescopes for SN1987A-like events, taking «det =100 
and different explosion energies. The shock breakout rate is 
calculated from the total core collapse SN rate per unit vol- 
ume, - 3 x 10" 4 Mpc" 3 yr" 1 dDahlenet al J 120121) . and the 
SN1987A-like eve nts fraction of that rate, which is ~ 3% 
JSmartt et al.ll2009h . 

Although SNellp comprise about 59% of the core collapse 
SN rate ( Smartt et al. 2009), their detection rate is not re- 
ported here as their breakout emission is expected to be ab- 
sorbed by the interstellar medium. In a similar fashion, SNelb 
comprise 10% of the core collapse SN rate, but the number of 
photons emitted at breakout is too low to detect. If the pro- 
genitors of this SN type are smaller than what was assumed 



10 



here, or the explosion energy is higher, than the emitted burst 
would peak at higher energies and might be detected by the 
BAT on-board Swift. 

As can be seen from table |3] XMM-Newton and ROSAT are 
the best candidates for detecting SN shock breakouts from 
SN1987A-like events, with > 30 events over the lifetime of 
each telescope. This means that there are expected to be many 
shock brea kout events in the archival d ata of each of these 
telescopes (Calzavara & Matzner (2004) co me to a similar 
conclu sion based on breakout calculations by Bli nnikov et al.l 
(2000)). We note however that previous transient surveys con- 
ducted with the archiva l data of ROSAT did not detect such 
bursts (IVikhlininlfT998t iGreineret ail 120001) . This poses a 
problem for the stellar envelope breakout scenario, as the pre- 
dictions of total emitted energy and total emitted photon num- 
ber are a robust result of the model, and mainly depend on the 
typical ejecta velocity and the progenitor's radius. Further- 
more, both the bolometric and temperature results are not sen- 
sitive to the density structure. Lack of shock breakout detec- 
tions may indicate that either the typical parameters assumed 
for the progenitors are grossly incorrect, or that the progen- 
itors are not characterized by a hydrostatic stellar envelope 
at the end of their lives (with a density profile monotonically 
decreasing towards the surface). 

6. PREVIOUS WORKS AND OBSERVATIONS 

We provide here a comparison to previous works in the lit- 
erature describing the temperature expected in shock break- 
outs from stellar envelopes, as well as a comparison to the 
proposed shock breakout candidate XRF080109. 

6.1. Temperature calculations in the literature 

Numerical works in the past have provided solutions for 
the emission following shock breakout from specific progen- 
itor profiles, but mostly for explosions with low ejecta ve- 
locity , A) < 0.05 (Ensman & Burrows 1992; IBlinnikov etaP 
2000). For example, Blinnikov et al. (2000) have calculated 
the expected breakout temperature for models describing the 
SN1987A progenitor, using pre-supernova density profiles 
which are based on stellar evolution models. However, these 
calculations did not incorporate incoherent Compton scatter- 
ing. 

We compare here our results with the ones reported for a 
BSG progenitor model with {M = l 4.67M e ,R = 4S . 5R Q } . 
For E = (0.72, 1 .02, 1 .34) x 10 51 erg, IBlinnikov et all d2000l) 
report a maximum bolometric luminosity of L max = 

(4.2.6.8.9.5) x 10 44 erg s" 1 , a maximum ejecta velocity of 
v m -dx = (2.8,3.3,3.7) x 10 9 cm s" 1 and a maximum tempera- 
ture of T « (90, 105, 1 15) eV, which roughly scales as L max . 
Our theoretical model for shock breakout from a gas with a 
power-law density profile and constat opacity provides sim- 
ilar results for the bolometric properties, but different ones 
for the temperature. Using the same set of parameters and 
taking v* = (E/M) l / 2 ,n = 0.34 and f p = 0.02, the maxi- 
mum bolometric luminosity and maximum ejecta velocity are 

(4.3.6.8.9.6) x 10 44 erg s" 1 and (2.5,3.1,3.6) x 10 9 cm s" 1 , 
respectively (see table 2 in Paper I and eqs. 7 & 25 in Pa- 
per II). Contrary to the bolometric properties, the difference 
in temperature is significant: according to eq. (l24l . the maxi- 
mum temperature is predicted to be T « (153,202,254) eV. 

Different predictions are obtained from the analysis of 
Nakar & Sari (2010), which provides analytical estimates for 



the breakout temperature, taking into account Comptoniza- 
tion of the bremsstrahlung emitted photons. However, their 
calculation assumes only local photon production, without ac- 
counting for diffusion or expansion. Using eq. 38 therein and 
for the same sets of parameters, the radiation is (marginally) 
out of thermal equilibrium, and the peak temperature is T « 
(95,175,274) eV. Therefore, our results are closer to the 
ones given by Naka r & Saril ([20 10), but the spread in max- 
imum tem perature between diff e rent sets of parameters ac- 
cording to iNakar & Saril (1201 Oh (IBlinnikov et al.l (120001) 1 is 
higher (lower) than our results which include incoherent scat- 
tering, photon diffusion and expansion. 

6.2. SN2008D as a shock breakout candidate 

X-ray flash (XRF) 080109/SN2 008D has been pro- 
posed as a shock breakout candidate (Soderberg et al. 2008; 
Chevalier & Fransson 2008; Katz et al. 2010), mainly based 
on its energetics (total released energy, typical luminosity). 
This XRF is characterized by a peak X-ray luminosity of 
10 43 erg s _1 in the 0.3-10 keV band and ~ 300 s duration, 
and a total emitted energy of E ~ 2.5 x 10 45 erg. 

The observed XRF fluence spectrum can be fitted by the 
universal fluence function (see eq. ( |29l and figure |5), with 
^peak = 3.9 keV, and is consistent with the breakout fluence 
spectrum up to photon energy of ~ 8 keV. The fit is somewhat 
in conflict with detection of photons above this energy, albeit 
with a small number of counts. In addition, an almost flat 
photon energy spectrum, S v ~ v° was suggested to agree with 
the observations ( Soderber g et al.| [2008; Modjaz et al. 2009), 
though it also has some tension with photons detected at ener- 
gies > 6 keV. If the emitted spectrum is indeed flat, then this 
event is not consistent with a shock breakout from a stellar 
envelope interpretation. Even without assuming a Wien spec- 
trum for the emission, the expected low energy tail of the flu- 
ence should be characterized by a positive spectral slope, as 
both the lu mino sity and the radiation temperature drop with 
time (see § 14. 11 1. A negative or flat slope would indicate that 
the emitted photons are not characterized by a specific radia- 
tion temperature, or that the main emission mechanism is ac- 
companied by a high energy power-law tail, both explanations 
are in disagreement with the stellar envelope shock breakout 
scenario. 

Assuming that the observed spectrum may be interpreted as 
a breakout spectrum, ignoring the above mentioned tensions 
between the model and the observations, i/ pea k = 3.9 keV im- 
plies, using eqs. d3~Tl i and d24l i. /3o ~ 0.12-0.17 (the varia- 
tion is due to the uncertainty in n p .o). This implies a terminal 
ejecta velocit y of v e j = (0.24—0. 34)c , which agrees with radio 
observations (Soderberg et al. 2008). Assuming k = 0.2, the 
progenitor's radius can then be found using the total emitted 
energy (see eq. 40 in Paper II) 



R-. 



1 kE 

hi^c 1 



1/2 



;5x 10 n cm, 



(44) 



where the correction of the emitted energy due to absorption 
is included and is only ~ 50%. The inferred radius is larger 
than a WR radius associated with this type lb SN, and is also 
a factor 5 larger than the value i nferred from UV/optical ob - 
servations in the spherical stage (Rabinak & Waxman 201 1). 
This radius also disagrees with the observed duration of the 
XRF, being a factor 20 smaller than the radius related to the 
duration, Ra = cAt w 10 13 cm. 



11 



If this event was indeed a shock breakout, the latter discrep- 
ancy can be explained either by assuming an non-spherical 
shock wave reaching different points on the surface at differ- 
ent times, or by assuming the stellar envelope is enshrouded 
by a moderately optically thick circumstellar material, with 
Tcsm *C (3q ' . The presence of a moderately optically thick ma- 
terial can change the duration of the observed burst, as this 
is determined by the light travel time across the radius of the 
last scattering surface. However, it does not significantly af- 
fect the total emitted energy or the fluence spectrum, as the 
extra mass is negligible, M c prn <C (k0o)~ 1 , and n o extra energy 
is added by the interaction (INakar & Sari 20l3). 



7. CONCLUSIONS 

Solutions for the temperature and photon number density 
were calculated for the problem of an RMS breaking out of 
a gas with a power law density profile, assuming constant 
Thomson scattering opacity, local Compton equilibrium and 
bremsstrahlung as the dominant emission rate. The radia- 
tion temperature is determined by three dimensional param- 
eters: breakout shock velocity vo, (pre-shock) density po and 
opacity k. Calculations presented here are the same for en- 
velopes with any hydrogen-helium composition, as the results 
depend on the atomic fractions of the elements %i and their 
atomic numbers Z, and weights A,- only through the combina- 
tion Ex t Zf fExiAi. 

The temperature at the surface, at an optical depth tvq/c = 
0, is similar to the temperature at an optical depth r = c/vq 
(see figure |2). The surface temperature rises with time to a 
sharp peak of width ~ to = c/(poVqK), at a time (0.1 - l)fo 
prior to maximum bolometric luminosity. The relation be- 
tween peak surface temperature, breakout shock velocity and 
density is well approximated by eq. (l24t . Tables[TJand[2]tabu- 
late the surface temperature as a function of time for different 
values of the parameter \ = ( v o/c)(po/(m p )/10 15 cm" 3 )" 1 / 36 . 
These tables can be used in order to calculate the temperature 
for different vo and po values not provided here, using linear 
interpolation in the parameter \ (see eq. d27l>). 

Assuming a Wien spectrum of the escaping photons, the 
instantaneous emission spectrum can be calculated using the 
bolometric luminosity and the surface temperature (see eq. 
d28ll). Accordingly, the luminosity in different energy bands 
can be calculated using eq. (132) . The spectral fluence, the 
integral over time of the spectral flux, shows a universal be- 
havior (see figure |5}, up to scaling of the fluence by £q = 
K~ l voc and of the photon energy by the peak surface temper- 
ature 7^eak- The predicted breakout fluence spectrum is non- 
thermal, characterized by a positive spectral slope in the low 
energy tail, and a peak energy at f pea k ~ 37^ e ak- 

The light curve in different energy bands for typical pro- 
genitors, corresponding to suggested progenitors of different 
SN types, was presented in figures |6]|8] It was shown that the 
prompt emission (on timescale R/c) following breakout from 
BSGs is mostly in the XRT band (0.3-10 keV), while for 
WRs it starts in the BAT band (15-150 keV) and then falls 
into the XRT band at times t k, R/(Avq). RSGs have a peak 



temperature lower than the energy band, but radiate signifi- 
cantly also in the XRT band (L ~ 10 44 erg s" 1 ). Our results 
differ from those of Nakar & Sari (2010), with a growing dis- 
crepancy as the radiation is driven further away from equilib- 
rium. In particular, the results of Nakar & Sari (2010) indicate 
that it would be difficult to observe RSG & WR breakouts in 
the XRT band, as their predicted temperatures for RSG and 
WR progenitors with parameters similar to those used here 
are closer to the UV and BAT bands respectively. 

We have compared our results to the numerical simulations 
of Blinnikov et al. (2000), and to the corresponding analytical 
estimates of Nakar & Sarii (120101) . for a given pre-supernova 
density profile assumed to describe the 1987 A progenitor. 
The peak bolometric luminosities and maximum ejecta ve- 
locities obtained in the numerical simulations are consistent 
with our results. The peak temperature in the numerical sim- 
ulations, which do not include incoherent Compton scatter- 
ing, roughly scales with maximum luminosity as Lj^, while 
the peak temperature in the analytical work of INakar & Saril 
(2010), which takes into account the Comptonization of pho- 
tons but not photon diffusion or the effect of hydrodynamic 
expansion, offers a bigger spread as a function of maximum 
luminosity. The results of this work are closer to the analytic 
estimates, but provide a smaller (larger) spread in the peak 
temperature as a function of maximum luminosity than the 
analytical (numerical) results. 

The detection rate of SN breakouts was evaluated for 
SN1987A-like events, considering the characteristics of dif- 
ferent X-ray telescopes (see table [3]calculated using eq. (l43l). 
> 30 SN1987A-like breakout detections are expected over the 
lifetimes of ROSAT and XMM-Newton. If such events are in- 
deed not detected it would indicate that the typical parame- 
ters assumed for the BSG progenitors are grossly incorrect, or 
that the progenitors' density profile does not monotonically 
decrease towards the surface, as expected for a hydrostatic 
stellar envelope. 

An analysis of the SN2008D/XRF 080109 event was per- 
formed, indicating that the observed fluence spectrum may 
be consistent with the breakout interpretation, provided that 
photons detected at energies > 8 keV are a statistical coin- 
cidence. Inferring the shock velocity from the peak fluence 
energy according to the universal breakout fluence, implies a 
terminal ejecta veloci ty which agrees with radio observations 
(Soderberg et al. 2008). On the other hand, the inferred radius 
of the event is R rs 5 x 10 11 cm, a value much higher than the 
expected radii of SNelb progenitors and also higher than the 
value inferred from optical /UV o bservations in the spherical 
stage (Rabinak & Waxman 2011), and much lower than the 
value suggested by the XRF duration. 

We thank David Burrows, Maryam Modjaz and Alicia 
Soderberg for useful discussions. N.S. and E.W are partially 
supported by UPBC and GIF grants. B.K. is supported by 
NASA through Einstein Postdoctoral Fellowship awarded by 
the Chandra X-ray Center, which is operated by the Smith- 
sonian Astrophysical Observatory for NASA under contract 
NAS8-03060. 



APPENDIX 

Numerical values ofT(t) 

Numerical results for the surface temperature as a function of time, T(t), for hydrogen-helium envelopes with n = 3 and 
n = 3/2 and several values of breakout shock velocity and pre-shock ion number density of n o = 10 15 cm" 3 are presented in 
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tables Q] and [2] The values presented are the surface temperature normalized by the peak surface temperature (which can be 
fitted by eq. (|24]|). and the time shifted by the time of peak surface temperature, normalized by to. In order to find the time 
dependent temperature for different values of shock velocity and density, the results can be linearly interpolated in the parameter 
X = /3o(An p ,o/10 15 cm" 3 )" 1 ' 35 , using eq. ( l27b . The time dependent temperature can then be used to determine the observed spectral 
luminosity (eq. d28Ti). fluence spectrum (eq. d29l i) and luminosity in different energy bands (eq. (l32l). 

Shock breakout detections by X-ray telescopes 

The expected number of detected SN shock breakout events from BSGs compatible with the progenitor of SN1987A, A^snsb, 
is presented for a number of active and non-active X-ray telescopes in table [3] For this analysis, the number of emitted photons 
for the different progenitors is taken from eq. (HTi and the number of detections is calculated using eq. (|43| >. assuming «,jet =100 
and taking the characteristics of the different telescopes presented in the table. 
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TABLE 1 

Normalized surface temperature for n ■- 
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TABLE 2 
Normalized surface temperature for n = 3/2 (table headings are the same as in tableQ} 
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"effective area for 1 keV photons 

b Field of view 

c Point spread function angular size, containing > 50% of the signal of a point source 

d Telescope's lifetime 

e Number of SN1987A-like shock breakout detections in the telescope's lifetime, according to eq. 1431 . with E m = 0.5 X 10 51 erg. 

= 10 51 erg. 

: 2 X 10 51 erg. 



f Number of SN1987A-like shock breakout detections in the telescope's lifetime, according to eq. 143 1 , with £;„ 
g Number of SN1987A-like shock breakout detections in the telescope's lifetime, according to eq. 1431 . with E m 



TABLE 3 
SN SHOCK BREAKOUT DETECTION WITH X-RAY TELESCOPES 



